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ABSTRACT 


Preliminary  studies  of  the  long-period  seismic  noise  and  the  atmospheric 
pressure  fields  were  made  at  the  Kleer  mine  near  Grand  Saline,  Texas,  during 
the  winter  of  1968.  The  results  of  the  studies  indicate  that  during  periods 
of  atmospheric  turbulence  the  vertical  and  horizontal  components  of  the 
seismic  noise  field  observed  at  the  surface  may  increase  in  total  power  by 
as  much  as  16  and  34  dB  within  the  period  range  from  20  to  60  seconds.  In 
contrast,  long-period  systems  operated  in  the  mine  show  little  change  in 
power  levels  regardless  of  the  atmospheric  conditions  at  the  suiface.  The 
observations  are  consistent  with  the  hypothesis  that  the  increase  observed 
at  the  surface  is  caused  by  atmospherically-generated  deformations  of  the 
earth.  However,  because  of  questions  regarding  the  surface  vaults  capability 
to  attenuate  pressure  changes,  the  possibility  that  the  noise  was  caused  by 
vault  leakage  must  be  considered  an  equally  likely  explanation  at  the  present 
time.  The  results  indicate  that  atmospherically-generated  seismic  noise  can 
seriously  degrade  the  capability  of  a  long-period  surface  installation  to 
detect  surface  waves.  Because  of  this  finding,  further  more  definitive  studies 
are  recommended. 


P REF ACt 


The  AR  method  and  its  derivatives  provide  a  reliable  discriminator  between 
seismic  signals  generated  by  large  magnitude  underground  nuclear  explosions 
and  earthquakes.  Unfortunately,  the  technique  cannot  be  applied  to  lower 
magnitude  events  because  of  a  relatively  high  surface  wave  detection  threshold. 
Since  the  identification  of  small  magnitude  events  is  cf  particular  importance 
to  undergound  nuclear  test  detection,  there  is  obviously  a  need  to  develop 
techniques  to  reduce  the  surface  wave  detection  threshold  or  equivalently  to 
attenuate  the  long-period  seismic  noise  field.  To  date  our  studies  have  been 
devoted  to  an  investigation  of  the  relationship  between  long-period  seismic 
noise  and  local  variations  in  the  atmospheric  pressure  field.  The  results 
of  these  studies  are  contained  in  this  report. 

This  report  is  a  compendium  of  three  related  studies.  The  first  study  describes 
the  observations  of  the  long-period  seismic  and  atmospheric  pressure  fields 
made  at  the  Kleer  Salt  Mine  near  Grand  Saline,  Texas.  The  basic  result  of 
this  study  is  that  the  long-period  noise  power  observed  at  the  surface  increases 
significantly  during  periods  of  atmospheric  turbulence. 

The  second  study  is  a  preliminary  theoretical  investigation  of  one  possible 
cause  for  the  observed  noise  increase.  In  this  paper,  the  response  of  an 
elastic  half  space  to  a  plane  pressure  wave  is  evaluated  in  order  to  determine 
the  conditions  under  which  the  pressure-generated  displacements  and  tilts  can 
contribute  significantly  to  the  seismic  background.  The  basic  result  of  this 
study  is  that  at  periods  greater  than  20  seconds  and,  in  media  with  moderately 
low  rigidities,  displacements  and  tilts  associated  with  pressure  waves  with 
amplitudes  on  the  order  of  100  microbars  can  contribute  significantly  to  the 
seismic  background.  The  results  further  indicated  that  if  the  speed  of  the 
pressure  wave  was  small  compared  to  the  speed  of  the  seismic  waves  in  the 
medium  the  deformation  could  be  described  by  the  mathematics  of  static  rather 
than  dynamic  elasticity. 

The  third  study  is  devoted  to  the  development  of  the  mathematics  necessary  to 
describe  the  static  deformation  of  a  layered  half  space  caused  by  an  arbitrarily 
distributed  surface  load. 

It  is  to  be  emphasized  that  the  results  presented  in  this  report  are  of  a 
preliminary  nature.  Observations  with  more  sophisticated  systems  are  presently 
being  made,  and  more  definitive  studies  of  the  relationship  between  long-period 
seismic  noise  and  local  atmospheric  pressure  fluctuations  are  underway.  The 
results  of  these  studies  will  be  published  in  a  later  report. 
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LONG-PERIOD  SEISMIC  NOISE  AND  ATMOSPHERIC  PRESSURE  VARIATIONS 


1.  INTRODUCTION 


Recent  theoretical  studies  by  Sorrells  (1969)  and  Der  (1969)  have  indicated 
that  variations  in  the  atmospheric  pressure  field  as  small  as  50-1U0  microbars 
can  generate  significant  seismic  noise  at  periods  greater  than  20  seconds. 
Experimental  observations  by  Capon  (1969)  at  the  Montana  LASA  tend  to  confirm 
the  theoretical  studies  and  imply  that  elastic  distortions  of  the  surface  of 
the  earth  induced  by  variations  in  atmospheric  pressure  are  indeed  a  major 
source  of  long-period  seismic  noise. 

During  the  winter  of  1968,  a  long-period  micorbarograph  and  two  long-period 
seismograph  systems  were  operated  at  the  Kleer  Mine  near  Grand  Saline,  Texas. 
The  records  produced  by  these  systems  were  examined  in  order  to  obtain  more 
information  regarding  the  possible  relationship  between  long-period  seismic 
noise  and  local  fluctuations  in  the  atmospheric  pressure  field.  The  purpose 
of  this  particular  report  is  to  document  the  results  of  the  preliminary  study, 
to  evaluate  their  significance,  and  to  outline  plans  for  further  study. 


2.  DATA  COLLECTION 


A  three- component  long- period  seismograph  system,  a  long-period  microbarograph, 
and  an  anemometer  were  operated  at  the  surface.  In  the  mine,  the  basic  sensor 
systems  included  not  only  a  three- component  long-period  system,  but  a  vertical 
and  horizontal  short-period  system  as  well.  Block  diagrams  of  the  surface 
and  mine  data  sensing  and  recording  systems  are  shown  in  figure  1.  The  system 
responses  of  the  long-period  seismic  and  microbarographic  systems  are  shown  in 
figures  2  and  3,  respectively.  The  seism! '  sensors  used  at  the  surface  were 
installed  in  typical  portable  system  vaults,  which  consist  of  a  hole  approxi¬ 
mately  2  feet  deep,  plywood-earth  retaining  walls  set  in  approximately  4  inches 
of  concrete.  The  space  around  the  seismometer  was  filled  with  fiberglass 
batting.  Then  the  vault  was  closed  with  a  plywood  cover,  two  layers  of  fiber¬ 
glass  batting,  and  a  sheet  of  clear  polyfilm.  The  seismometers  used  in  the 
mine  were  set  on  the  floor.  Both  the  mine  and  surface  long-period  vertical 
seismometer  cases  were  sealed  while  the  horizontal  seismometer  cases  were  not. 
The  vertical  separation  of  the  mine  and  surface  systems  was  approximately  600 
feet  while  as  shown  in  figure  4  the  horizontal  separation  was  approximately 
400  feet. 
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Figure  1.  Block  diagram  of  the  instrumentation  systems  operated  at  Grand  Saline,  Texas 
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Figure  3 


System  response  of  the  long-period  microbarograph 
operated  at  Grand  Saline,  Texas 
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Figure  4.  Sketch  map  of  the  Kleer  mine  showing  the  relative 
location  of  the  mine  and  surface  systems 
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3.  DATA  ANALYSIS 


Power  density  spectra  for  individual  records  and  coherency  spectra  for 
selected  pairs  of  records  were  the  basic  analytical  tools  used  during  this 
investigation.  These  were  computed  from  record  segments  chosen  to  represent 
various  conditions  of  atmospheric  pressure  variations.  A  total  of  six  such 
segments,  each  approximately  30  minutes  in  duration,  were  selected  for  the 
preliminary  analysis.  The  selection  was  made  from  16-mm  film  reproduction  of 
the  data  stored  on  magnetic  tape.  Within  each  of  the  selected  time  intervals, 
the  outputs  of  both  sets  of  seismographs  and  the  microbarograph  were  low-pass 
filtered  to  prevent  aliasing,  digitized  at  the  rate  of  two  samples  per  second, 
then  multiplexed  and  stored  on  digital  tapes.  Prior  to  computing  the  power 
and  coherency  spectra,  the  seismic  data  were  normalized  by  dividing  by  the 
magnification  at  25  seconds;  and  linear  trends  were  removed  from  both  the 
seismic  and  microbarographic  records.  Both  the  power  density  and  cross  power 
density  spectra  were  computed  from  blocks  of  3700  points  using  127  lags  and 
Parzen  smoothing.  Power  density  spectra  were  computed  for  all  records.  Cross 
power  density,  coherency  and  phase  spectra  were  computed  for  corresponding 
seismic  records  taken  on  the  surface  and  in  the  mine.  In  addition,  cross 
power  density,  coherency,  and  phase  spectra  were  computed  for  records  taken 
from  the  microbarograph  and  the  surface  horizontal  seismographs.  Examples 
of  the  results  of  the  computation  are  displayed  anl  discussed  in  the  next 
section . 


4.  DISCUSSION  OF  RESULTS 


4.1  RECORDS  TAKEN  DURING  CALM  PERIODS 

The  power  spectra  shown  in  figure  5  are  typical  of  the  seismic  data  recorded 
on  the  vertical  seismographs  during  calm  periods.  Note  that  for  periods  less 
than  about  20  seconds,  the  spectral  content  of  the  noise  observed  on  the 
surface  and  in  the  mine  is  virtually  the  same.  This  was  not  unexpected  since 
the  seismic  noise  in  the  period  range  from  6-20  seconds  is  known  to  consist 
primarily  of  fundamental  mode  Rayleigh  waves  with  wavelengths  much  greater 
than  the  distance  separating  the  two  observation  points.  The  coherency  measure¬ 
ments  shown  in  figure  6  confirm  that  within  this  period  range  the  noise  recorded 
by  both  systems  is  identical.  At  the  longer  periods,  however,  the  situation  is 
quite  different.  As  seen  in  figure  6,  there  is  a  progressive  reduction  in  the 
coherence  between  the  noise  observed  on  the  surface  and  within  the  mine.  This 
is  accompanied  by  a  monotonic  increase  in  the  noise  power  observed  on  the  sur¬ 
face  relative  to  that  observed  in  the  mine  (see  figure  7) .  Clearly,  there  is 
no  apparent  relationship  between  the  seismic  noise  observed  by  the  two  systems 
within  this  period  range. 
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Figure  5 


Power  spectra  typical  of  the  data  recorded  in  the  mine 
and  on  the  surface  by  the  vertical  seismographs  during 
a  calm  period.  Spectra  are  not  corrected  for  instrument 
response . 
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The  data  recorded  by  the  horizontal  seismographs  during  calm  periods  show  some 
characteristics  similar  to  those  described  in  the  previous  paragraph.  An 
example  of  the  power  spectra  computed  from  data  taken  from  the  transverse 
horizontal  seismographs  located  on  the  surface  and  in  the  mine  are  shown  in 
figure  8.  These  were  computed  from  data  taken  during  the  same  interval  of 
time  as  the  vertical  data  discussed  previously.  The  coherence  between  the  two 
sets  of  horizontal  data  is  shown  in  figure  9.  Like  the  vertical  data  the 
horizontal  spectra  may  be  divided  into  two  period  ranges;  a  short-period  range 
characterized  by  high  coherence  in  which  the  power  levels  observed  upon  the 
surface  and  within  the  mine  are  the  same,  and  a  long-period  range  within  which 
the  two  data  sets  are  incoherent  and  the  noise  power  observed  upon  the  surface 
is  considerably  higher  than  that  observed  in  the  mine.  The  principal  difference 
is  that,  whereas,  the  15-second  microseism  peak  is  clearly  visible  on  both 
vertical  spectra,  it  is  almost  completely  obscured  on  the  horizontal  spectra. 

The  interfering  fields  as  recorded  by  the  two  horizontal  seismographs  are 
apparently  generated  by  independent  sources  as  shown  by  the  reduction  in  the 
coherence  of  the  15-second  peak.  The  fact  that  even  during  calm  periods  the 
noise  at  periods  greater  than  about  15  to  20  seconds  is  considerably  higher 
on  the  surface  than  in  the  mine  merits  further  consideration.  It  cannot  be 
related  to  variable  atmospheric  conditions  because  the  samples  analyzed  were 
specifically  selected  from  intervals  of  time  during  which  there  were  no  measure- 
able  wind  and  microbaric  activity  was  minimal.  This  conclusion  is  further 
reinforced  by  the  fact  that  there  is  no  coherence  between  data  recorded  by  the 
microuarograph  and  by  the  vertical  and  horizontal  seismographs  (see  figure  10) . 
There  are  several  possible  explanations  for  the  existence!  of  the  high  noise 
at  the  surface.  First  of  all,  the  noise  could  be  system  generated.  This  seems 
unlikely  since  it  is  found  at  a  relatively  high  level  only  on  the  surface 
systems.  A  more  likely  explanation  is  that  it  is  in  some  way  associated  with 
the  type  of  installations  used  on  the  surface.  It  will  be  recalled  from 
section  2  that  while  the  surface  vaults  were  covered,  no  specific  attempts  were 
made  to  tightly  seal  them.  Thus,  it  is  possible  that  some  of  the  noise  observed 
at  the  longer  periods  could  be  related  to  disturbances  within  the  vault. 

Another  possibility  is  that  the  noise  is  related  to  local  surface  activity 
associated  with  the  operation  of  the  mine.  No  conscious  effort  was  made 
during  this  particular  investigation  to  select  samples  during  periods  when 
activity  associated  with  the  mine  was  likely  to  be  low.  Thus,  it  is  possible 
that  the  samples  could  be  contaminated  to  some  decree  by  noise  generated  through 
local  surface  activity.  This  possibility  will  be  thoroughly  examined  during 
a  forthcoming  investigation. 


4.2  RECORDS  TAKEN  DURING  TURBULENT  PERIODS 

The  power  spectra  shown  in  figure  11  and  figure  12  display  features  that  are 
typical  of  the  data  recorded  by  the  vertical  and  horizontal  seismographs  during 
turbulent  periods.  The  basic  difference  between  these  spectra  and  those  taken 
during  calm  periods  is  the  relatively  large  increase  in  the  noise  power  observed 
at  the  surface.  The  nature  of  this  increase  may  be  seen  more  clearly  from 
figure  13.  In  this  figure,  the  ratio  between  the  spectra  observed  at  the 
surface  during  a  calm  and  turbulent  interval  has  been  plotted  as  a  function  of 
period  in  the  range  from  20  to  64  seconds.  Note  that  both  the  surface  vertical 
and  transverse  ratios  show  a  trend  similar  to  the  microbarographic  ratio.  All 
ratios  tend  to  decrease  in  a  regular  manner  as  a  function  of  increasing  period. 
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Figure  8.  Power  spectra  typical  of  data  recorded  in  the  mine  and 
on  the  surface  by  the  transverse  horizontal  seismographs 
during  a  calm  period.  Spectra  are  not  corrected  for  system 
response . 
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Figure  9.  The  coherence  between  data  recorded  by  the  transverse  horizontal  seismographs 
located  in  the  mine  and  on  the  surface  -  data  recorded  during  a  calm  period 
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located  at  the  surface  and  the  microbarograph  -  data  recorded  during 
calm  interval 


POWER  DENSITY 


Figure  11.  Power  spectra  computed  from  data  taken  from  mine  and 

surface  vertical  seismographs  during  a  turbulent  period. 
Spectra  are  not  corrected  for  system  response.  G  5821 
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Figure  13.  The  ratio  of  the  power  spectra  computed  from  data  observed 
at  the  surface  during  a  turbulent  period  relative  to  those 
computed  from  data  recorded  during  a  calm  period 
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However,  both  the  values  of  the  microbarographic  ratios  and  their  rate  of 
decrease  with  increasing  period  is  greater  than  for  the  seismic  data.  The 
results  shown  in  this  figure  strongly  suggest  that  the  increase  in  the  long- 
period  noise  field  observed  at  the  surface  during  turbulent  intervals  is 
directly  related  to  increases  in  the  atmospheric  pressure  field.  Additional 
support  for  this  contention  comes  from  estimates  of  the  degree  of  coherence 
between  the  data  recorded  by  the  microbarograph  and  the  surface  seismographs. 
These  estimates  are  plotted  as  a  function  of  increasing  period  in  figure  14. 
Notice  that  at  periods  greater  than  about  20  seconds,  a  significant  fraction 
of  the  total  noise  power  observed  by  the  transverse  seismograph  can  be 
correlated  atmospheric  pressure  variations.  The  degree  of  coherence  between 
the  microbarograph  and  the  surface  vertical  is  considerably  smaller,  indicating 
that  the  percentage  of  atmospherically-generated  noise  relative  to  the  total 
field  observed  by  this  system  is  much  lower  than  by  that  observed  by  the 
surface  transverse.  This  is  consistent  with  the  observations  given  in 
figure  13,  which  show  that  the  increase  in  the  power  of  the  vertical  component 

of  the  noise  field  is  about  6  dB  less  than  the  increase  in  the  horizontal 

component . 

In  contrast  to  the  seismic  data  observed  at  the  surface,  there  is  little 
significant  change  in  the  data  recorded  in  the  mine  regardless  of  the  atmo¬ 
spheric  conditions.  This  may  be  seen  from  figure  15,  in  which  the  ratio  of  the 

power  observed  in  the  mine  during  calm  and  turbulent  periods  is  plotted  as  a 
function  of  increasing  period.  The  shaded  area  shown  in  this  figure  indicates 
the  90  percent  confidence  interval.  Ratios  which  fall  in  this  range  cannot 
be  considered  significantly  different  from  1  (0  dB)  .  Note  that  there  is  no 
significant  difference  in  the  power  levels  observed  by  the  transverse  hori¬ 
zontal  seismographs  in  the  period  range  from  20  to  64  seconds.  The  vertical 
component  of  the  seismic  noise  field  observed  in  the  mine  shows  a  small  in¬ 
crease  at  periods  less  than  about  30  seconds,  but  the  increase  is  about  10  dB 
less  than  that  observed  at  the  surface.  It  is  clear  from  these  results  that 
atmospheric  turbulence  at  the  surface  does  not  contribute  significantly  to  the 
seismic  noise  field  in  the  mine. 

In  addition  to  the  specific  features  described  in  the  previous  paragraphs, 
there  is  a  general  tendency  for  the  seismic  noise  observed  at  the  surface  to 
increase  as  the  noise  recorded  by  the  microbarograph  increases.  This  tendency 
is  illustrated  in  figures  16a  and  16b.  The  diagrams  were  constructed  in  the 
following  manner.  First,  both  the  sjismic  and  the  microbarographic  spectra 
collected  during  various  intervals  of  atmospheric  turbulence  were  normalized 
to  a  common  base  by  dividing  them  by  the  spectra  collected  during  a  calm  period. 
Then,  the  normalized  spectra  were  averaged  over  the  period  range  extending 
from  21.3  to  64  seconds.  The  resulting  number  is  a  measure  of  the  relative 
mean  change  in  long-period  power,  which  pccurred  during  that  interval  as  com¬ 
pared  to  the  power  observed  during  a  calm  interval.  Then  in  order  to  delineate 
any  possible  relationship  between  changes  in  microbarographic  power  and  changes 
in  seismic  power,  the  mean  seismic  values  that  were  computed  from  data  collected 
at  the  surface  were  plotted  against  the  mean  microbarographic  power  change 
characteristic  of  that  interval.  As  shown  in  figures  16a  and  16b,  both  the 
ho7izontal  and  vertical  components  of  the  surface  seismic  field  tend  to  increase 
monotonically  as  the  microbarographic  power  increases,  but  not  at  the  same  rate. 
The  horizontal  rate  of  increase  appears  to  be  considerably  larger  than  the 
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Figure  14.  The  coherence  between  the  seismic  data  recorded  at  the  surface 

and  the  microbarograph  during  a  turbulent  period  g  5824 
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Figure  15.  The  ratio  of  the  power  spectra  computed  from  data  observed  in  the 

mine  during  a  turbulent  period  relative  to  those  computed  from  data 
recorded  during  a  calm  period 


G  5825 
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vert^ai.  If  similar  seismic  quantities  computed  from  the  data  collected  in 
the  mine  are  plotted  against  relative  increases  in  microbarographic  power,  the 
results  shown  in  figures  17a  and  17b  are  obtained.  These  daw  show  ?hat 
atmospheric  pressure  fluctuations  at  the  surface  have  no  apparent  effect  on 
the  horizontal  component  of  the  long-period  seismic  noise  field  observed  in 
the  mine  The  vertical  component  does  show  a  weak  tendency  to  increase  as 
the  microbarographic  power  increases;  however,  the  rate  of  increase  is  much 
less  than  that  observed  on  the  surface. 


5.  SUMMARY  OF  RESULTS 


The  basic  results  of  this 
following  manner. 


preliminary  investigation  may  be  summarized  in  the 


a.  During  calm  intervals  and  in  the  period  range  from  about  6  to  20 
seconds,  the  power  levels  recorded  in  the  mine  and  on  the  surface  are  quite 
similar  and  show  a  high  degree  of  coherence  indicating  identical  sources.  At 
periods  greater  than  20  seconds,  the  coherence  between  the  data  recorded  at 
the  surface  and  within  the  mine  decays  rapidly  and  the  power  levels  observed 
at  the  curface  increase  dramatically  relative  to  those  observed  within  the 

mine.  The  cause  of  the  higher  noise  at  the  surface  is  not  known  at  the  present 
time.  r 


b.  During  periods  of  atmospheric  turbulence,  the  seismic  noise  field 
observed  at  the  surface  contains  a  component  which  is  related  to  fluctuations 
in  the  atmospheric  pressure.  The  addition  of  this  component  to  the  surface 
seismic  field  results  in  a  general  tendency  for  the  total  field  to  increase 
as  the  atmospheric  pressure  fluctuations  intensify.  This  component  appears 

to  be  absent  or  at  least  strongly  attenuated  in  the  seismic  field  observed 
in  the  mine. 


6.  EXPLANATION  OF  RESULTS 


As  stated  previously,  the  cause  of  the  relatively  high  noise  levels  observed 
at  the  surface  during  calm  periods  is  not  presently  known.  It  could  be  caused 
by  a  variety  of  factors,  none  of  which  were  specifically  examined  during  the 
course  of  this  investigation.  Among  these  are  the  presence  f  higher  system 
noise  on  the  units  operated  at  the  surface;  unspecified  local  noise  associated 
with  the  nature  of  the  surface  installations;  and  noise  created  by  activities 
associated  with  the  operation  of  the  mine. 
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MICROBAROMETRIC  POWER  (dB) 


Figure  17.  The  mean  increase  in  seismic  noise  power  level  observed  in  the 
mine  versus  the  mean  increase  in  the  microbarographic  noise 
power  level 

G  3827 

-22- 


TR  70-12 


The  results  presented  display  several  attributes  which  suggest  that  during 
turbuient  periods  a  significant  fraction  of  the  surface  noise  field  could^e 

variations  of  the  surface  caused  by  atmospheric  pressure 

variations.  The  strongest  arguments  in  favor  of  this  interpretation  are: 

a.  There  is  a  significant  reduction  in  the  atmospherically-generated 
component  of  the  long-period  seismic  noise  field  observed  in  the  mine.  This  is 
in  qualitative  agreement  with  earlier  calculations  by  Sorrells  (1969)  which 

nnniC^!d  ^  **  ^  n°iSG  Were  caused  by  variable  atmospheric  loads  *  imposed 
p  he  surface  then  at  the  depth  at  which  the  seismic  systems  were  located 

(apprcnamateiy  600  feet  below  the  surface)  it  should  be  reduced  by  roughly 
40  dB  relative  to  its  strength  at  tne  surface.  Consequently,  the  general  lack 
of  an  atmospherically-related  component  is  consistent  with  the  hypothesis  of 
elastic  loading  by  atmospheric  pressure. 

b.  Durings  periods  of  increased  atmospheric  noise  power,  the  noise 
observed  by  the  horizontal  system  increases  more  than  the  noise  observed  by 
the  vertical  system.  This  observation  agrees  qualitatively  with  Sorrells' 
results,  which  indicate  that  at  periods  greater  than  about  20  seconds  a  given 
pressure  change  will  generate  apparent  horizontal  displacements  due  to  tilt 
which  are  at  least  an  order  of  magnitude  greater  than  the  corresponding  vertical 
displacements.  Thus,  if  the  increase  in  the  seismic  noise  field  observed  at 
the  surface  were  caused  by  pressure-generated  distortions  of  the  surface  then 
the  increase  in  noise  power  should  be  much  greater  on  the  outputs  of  the*hori- 
zontal  systems  than  on  the  vertical. 

It  must  be  stressed  that  these  observations  do  not  prove  the  hypothesis  that 
the  noise  increase  observed  during  turbulent  intervals  is  caused  by  pressure- 
generated  distortions  of  the  earth.  They  merely  argue  for  its  plausibility. 
Until  steps  can  be  taken  to  assure  that  the  surface  vaults  do  not  pass  atmo¬ 
spheric  pressure  variations  in  the  period  range  from  about  20-60  seconds  the 
possibility  that  the  noise  increase  is  caused  by  pressure  variations  within  the 
vaults  must  be  considered  an  equally  likely  explanation 


7.  SIGNIFICANCE  OF  RESULTS 


The  observations  reported  have  a  direct  bearing  on  the  problem  of  detecting 
long-period  surface  waves.  It  has  been  shown  that  during  intervals  of  moderate 
atmospheric  turbulence  and  within  the  period  range  from  20-60  seconds  the 
vertical  component  of  the  surface  seismic  noise  field  can  increase  by  as  much 
as  16  dB  in  power,  while  the  horizontal  component  can  increase  by  as  much  as 
34  dB  in  power.  These  numbers  imply  that  a  station's  capability  to  detect 
surface  waves  can  be  seriously  degraded  during  periods  of  atmospheric  turbulence. 
This  phenomena  is  by  no  means  isolated  to  the  Grand  Saline  area.  Capon  (1969) 
has  shown  that  at  the  Montana  LASA  and  within  the  ptriod  range  20-40  seconds 
atmospherically-generated  noise  makes  up  at  least  40  percent  of  the  total  power 
at  least  50  percent  of  the  time.  y 
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It  also  has  been  shown  <- that  the  noise  levels  observed  within  the  mine  stay 
essentially  the  same  regardless  of  the  atmospheric  conditions  at  the  surface. 

This  phenomena  is  apparently  related  to  the  fact  that  the  seismic  disturbances 
created  by  variable  atmospheric  loading  dissipate  very  rapidly  with  depth  as 
predicted  by  Sorrells  (1969).  Van  Saggern  (1970)  has  reported  similar  observa¬ 
tions  at  Murphy  Dome,  Alaska,  so  that  the  absence  of  atmospherically-generated 
seismic  noise  at  moderate  depths  is  not  unique  to  Grand  Saline. 

These  observations  when  taken  together  imply  that  the  surface  wave  detection 
capability  of  a  given  site  or  array  cf  sites  can  be  significantly  improved  by 
locating  the  sensors  at  moderate  depths .  This  improvement  would  be  the  result 
of  two  factors.  First  of  all,  the  observations  documented  in  this  report 
indicate  that  the  vertical  component  of  the  noise  field  would  be  reduced  by  as 
much  as  6-8  dB  in  amplitude  and  the  horizontal  components  would  be  reduced  by 
as  much  as  17  dB  in  amplitude.  Secondly,  after  removal  of  the  atmospherical lv- 
generated  fraction,  the  remainder  of  the  noise  field,  at  least  in  the  period 
range  from  20-40  seconds,  is  anticipated  to  be  coherent  over  relatively  large 
distances  (Capon,  1969).  Thus,  arrays  and  multichannel  filtering  techniques 
could  be  employed  to  achieve  a  further  reduction  of  about  6  dB.  Therefore, 
the  total  increase  in  surface  wave  signal-to-noise  ratio  could  be  as  much  as 
12-14  dB.  If  such  a  gain  could  be  realized  it  would  raise  the  surface  wave 
detection  capability  to  s  level  comparable  to  the  present  short-period  P  wave 
detection  capability. 

One  further  point  should  bt  made.  Recently  there  has  been  considerable  interest 
in  altering  the  response  of  long-period  systems  so  that  reliable  observations  of 
seismic  data  with  periods  up  to  100  seconds  can  be  made.  The  observations 
reported  in  the  previous  paragraphs  are  pertinent  to  this  problem  in  that  they 
indicate  that  a  significant  fraction  of  the  noise  in  the  period  range  extending 
from  20  seconds  is  related  to  atmospheric  turbulence.  It  is  anticipated  that 
this  noise  will  become  even  more  intense  as  the  period  range  is  extended.  The 
reasons  for  this  are  as  follows: 

a.  It  has  been  generally  observed  that  the  atmospheric  pressure  increases 
linearly  with  increasing  period. 

b.  The  earth's  vertical  displacement  response  to  atmospheric  pressure 
variations  also  increases  linearly  with  the  period  (Sorrells,  1969). 

c.  The  horizontal  seismometer's  tilt  response  increases  as  the  square 
of  the  period  (Rodgers,  1968). 

Since  significant  atmospherically-related  seismic  noise  has  already  been  observed 
at  the  shorter  periods,  it  can  be  concluded  from  the  statements  listed  above 
that  it  will  constitute  an  even  greater  fraction  of  the  noise  observed  at  the 
longer  period  ranges.  Thus,  it  appears  likely  that  in  order  to  achieve  high 
magnifications  with  the  extended  range  long-period  systems  they  must  be  operated 
at  depths  on  the  order  of  200  to  300  meters  or  greater. 
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8.  RECOMMENDATIONS 


In  view  of  the  significan 
of  surface  wave  detection 
Grand  Saline.  The  basic 


ce  of  these  preliminary  observations  to  the  problems 
*  it  is  recommended  that  the  studies  be  continued  at 
objectives  of  these  studios  will  be: 


a.  To  determine  the  source  of  the  relatively  high  noise 
at  the  surface  during  calm  periods  aid,  if  possible,  eliminate; 


levels 


observed 


b.  To  estimate  the  fraction  of 
extending  from  about  30  to  100  seconds 
atmospherically-generated  deformations 


seismic  noise  power,  within  a  pass  band 
,  which  can  be  attributed  directly  to 
of  the  earth. 


In  order  to  satisfy  these  basic  objectives, 
accomplished: 


the  following  tasks  must  be 


a  Seismic  systems  that  possess  relatively  high  magnifications  in  the 

100  f Conds  mu3t  be  designed,  fabricated,  tested,  and 
installed  both  on  the  surface  and  within  the  mine  at  Grand  Saline. 

i  _  b*  1116  vaults  containing  the  seismometers  must  be  sealed  to  the  extent 
that  pressure  variations  with  periods  as  great  as  400  seconds  are  strongly 
attenuated  within  the  vault. 


c. 


A  microbarograph  array  must  be  operated  at  the  surface  in  close 
proximity  to  the  seismometer  vaults  in  order  to  determine  the  frequency  wave 
number  content  of  the  atmospheric  noise. 


In  the  event  that  atmospherically-generated  seismic  noise  does  prove  to  contri¬ 
bute  significantly  to  the  surface  seismic  background  within  the  pass  band  of 
interest,  methods  to  reduce  its  effect  must  be  developed  and  evaluated  At 
the  present  time,  it  appears  that  this  contribution  could  be  significantly 
reduced  by  either  of  the  following  methods,  *  7 


a. 


Installation  of  the  sensors  at  moderate  depths; 


b.  Estimation  of  the  atmospherically-generated  component  of  the  seismic 
noise  using  microbarograph  array  data.  This  estimation  would  then  be  used  to 
eliminate  the  actual  component  by  subtracting  it  from  the  observed  field. 

Both  methods  should  be  evaluated  thoroughly.  It  is  recommended  that  the  estima¬ 
tion  process  be  evaluated  during  the  studies  at  Grand  Saline  since  high  quality 
microbarographic  data  will  be  available.  The  degree  of  attenuation  with  depth 
however,  cannot  be  studied  at  Grand  Saline  since  it  is  difficult  if  not 
impossible,  to  move  the  seismic  system  vertically  within  the  mine.  A  more 
effective  technique  for  evaluating  this  method  would  be  to  operate  a  long-period 
tnaxial  seismograph  system  at  various  depths  within  a  borehole  and  compare  the 
results  with  observations  taken  at  the  surface.  More  specific  recommendations 
regarding  a  study  that  will  employ  borehole  observations  to  study  the  attenuation 
of  atmospherically-generated  noise  with  depth  will  be  made  at  a  later  date  pendin 
the  outcome  of  the  studies  at  Grand  Saline.  P 
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THE  RESPONSE  OF  AN  ISOTROPIC  HALF  SPACE  TO 
_ _ A  PLANE  PRESSURE  WAVE 


1.  INTRODUCTION 


A  recent  article  by  Capon,  1968,  has  suggested  that  seismic  noise  created  by 
local  variations  in  the  atmospheric  load  may  contribute  significantly  to  the 
long-period  noise  field  at  periods  greater  than  20  seconds.  This  note  is  the 
first  in  a  series  of  short  papers  describing  the  results  of  theoretical  and 
experimental  investigations  of  this  hypothesis.  This  particular  report  is 
devoted  to  an  analysis  of  the  response  of  a  perfectly  elastic,  isotropic 
half  space  to  a  plane  pressure  wave  which  moves  at  a  constant  rate  well  below 
the  seismic  wave  speeds  of  the  medium.  This  is  a  relatively  simple  problem 
which  can  be  solved  without  resorting  to  numerical  techniques.  It  has  been 
considered  by  various  authors  (cf.  Cole  and  Heath,  1958;  Papadou^as,  1963; 
Ang,  1960;  Payton,  1967) .  Our  purpose  in  re-examining  the  solution  is  to 
determine  if  the  earth's  response  to  small  changes  in  atmospheric  pressure  is 
sufficiently  sensitive  to  warrant  the  study  of  more  sophisticated  models. 


2.  THEORY 


Let  x  and  z  denote  cartesian  coordinates  in  directions  parallel  and  normal 
to  the  surface  of  a  half  space.  Let  the  plane  z  =  0  coincide  with  the  surface 
of  the  half  space  and  take  z  positive  into  the  medium  and  x  positive  to  the 
right  of  an  arbitrarily  chosen  origin  on  the  surface.  Let  ct  and  8  denote 
the  compressional  and  shear  wave  velocities  of  the  medium  and  let  p  denote 
the  density.  Suppose  that  a  pressure  disturbance  p  (x,t)  is  imposed  upon  the 
surface  of  the  half  space.  Then  letting  azz(x,z,t)  and  azx  (x,z,t)  denote 
the  normal  and  tangential  stresses,  the  boundary  conditions  are 


°zz  °»  ^  =  P  Cx»  D 


u>) 


,i“t  Q-ikx 


dkdic 


CD 


°zx  Cx»  0,  t)  =  0  (2) 

Now  define  displacement  potentials  0  and  ¥  by 

9  (x,  z,  t)  =  jj  A  0,k)  eiut  e-kxe'YZ  dkdu  (3) 
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where  A  and  B  are  functions  to  be  determined  by  the  boundary  conditions.  It 
may  be  readily  verified  that  these  potentials  satisfy  the  wave  equations 

C5) 

a2  9  t2 

V2  *  =  1  a2  (6) 


provided  that 


Y 


2 


(7) 


where 


(8) 


(9) 


k 


2 


(10) 


In  terms  of  the  displacement  potentials  the  normal  and  tangential  stresses 
become  (Ewing  et  al,  1957,  p.  26) 


2  9 2  <p  +  d2f 

°zz  (x,z,t)  =  X  V  $  +  2y  (  9z2  9x9z  ) 

2  2  2 

3  <p  +  31  -  3  y 

°zx  =  ^(^  3x3z  3x^  3z2) 
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where 


X  =  p  (a2-2  02) 


(13) 


y  =  p 


(14) 


Then  using  the  boundary  conditions  defined  by  equations  1  and  2  in  equations  11 
and  12  together  with  the  definitions  given  by  3  and  4  the  following  relations  for 
A  and  B  are  obtained 


A  Go ,k)  =  isiidal  P  (15) 

y  F(a),k) 


B  (<o,k)  =  2i  k  Y  P(u>,k) 

y  F  (a),k)  (16) 

where 

F(o),k)=  (2k2  -  k2  )2  -  4k2  yy’  (17) 


Now  let  U(x,z,t),  W(x,z,t)  denote  the  displacements  in  the  half  space  in  the 
x  and  z  directions,  respectively.  In  terms  of  the  potentials 


u  =  3<j>  - 

3j> 

(18) 

3x 

3z 

W  =  34>  + 

3j> 

(19) 

3z 

3x 

and  using  equations  15  and  16  to  replace  A  and  B  in  equations  3  and  4  and  then 
performing  the  operations  described  by  equations  18  and  19,  the  following 
expressions  for  U  and  W  are  obtained 


U(x,z,t) 


P  (k,a>)  k  ,[(2k2-k8 
F(k,uO 

eiu)t  e-ikx  dkdu) 


)e"Yz-2YY’e~r  z] 


(20) 
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00  00 


W(x,z,t)  =  /  /  P  (k,0))  y  [(2k2-k62)e‘YZ-2  kVY’2i 

4n2uy  J  F(k,u>) 


(21) 


.00  .00 


eiwte-ikx  dkda) 


2.1  RESPONSE  TO  A  PLANE  WAVE  SOURCE 

The  remainder  of  this  note  will  be  devoted  to  a  study  of  the  half  space  response 
to  plane  wave  sources.  A  plane  wave  source  is  defined  by  the  relation 


P  (x,t)  =  p  (t  -_2L) 

c 


(22) 


where  c  is  the  speed  of  the  plane  wave  and, for  the  purposes  of  this  study,  is 
restricted  to  values  such  that  c  <<6.  Now 


00  00 

•«-// 


p(x,t)  e'iUt  eikx  dt  dx 


.00  .00 


Therefore  in  the  case  of  the  plane  wave  source 


;>u))  =  j  j  p(t-x)  e“iu)t  eikx  dt  dx 


.00  .00 


Let  £  =  t  -  21 
c 


00  00 

,„)  =  j  ei  x  j 


P  (£)  e"iu,5d£  dx 


(25) 


=  2ttP(uj)  6  (k-u)) 


(26) 


Substitution  of  equation  26  into  equations  20  and  21  and  integration  over 
k  yields 
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00 


U(x,z,t)  = 


pm  I  1(2  J  •^)e-rW!.2)Ml,M/W‘] 


F  (oj  ,ai  ) 


elu(t-x)  . 
e  c  d" 


W(x,z,t)  =  /  p(<*0Y  (w) 


[(2  \ 
c* 


£  •* , 


-Y'C“)Z 


where 


F  (  «,<»>) 


.i^Ct-x)  dw 


Y  (.<*>) 


C  2 

n 


2,1/2 


(29) 


JlLd-  cV/2 

y'M=  c 


(30) 


Using  the  absolute  value  of  w  assures  that  the  displacements  vanish  as  z  - 
Now  for  the  case  where  c  <<  8  we  may  write 


lid.  (1  -  -  °2)  +  Of-^H 

Y  O)  *  c  L  2  ^2J  UC  4' 


Therefore 


1  2 
1  c 


Y*  (“Osr  c  C1  *  J  “T)  +  or—) 

e2  V 


(32) 
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♦  J"L  ]e-Mz 


(36  cont'd) 


e2  2 
~  e~a 
*  “2T 
a  3 


a 


Similarly  part  of  equation  (28)  can  be  approximated 


Y  0*0 


ad) 

/ e 


.y 


e-y’Mz  j 


2 


c2 

7 


(37) 


e2-a2  M3 
a2  82  c 


and 


CO 


F(^.»a))  * 


4 

in _ 

c4 


[(2-  Si)2-  4  (1  1  Sl)(l  -  1  4  ] 

e2  “ 


8' 


[4  -  4_ci 

71 


-  4  (i  -  i  (£♦£))  ] 
«2  62 


(38) 


r  tdz* 

a2  8 


(39) 
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Now 


02-a2 


(40) 


a 


2 


X+2jj 
A  +  u 


(41) 


Using  equations  3b  through  41  in  27  and  28,  one  obtains  for  U  and  W 


U(x,z,t) 


00 


-00 


eia)  (t~c  )  dui 


(42) 


W  (x,z,t) 


P  (<*>)(• 


A+2y 

A+V 


M_z 

c  J  e" 


"  z 


,  iu(t  - 


*)  , 
c  Oo) 


0) 


(43) 


Equations  42  and  43  describe  a  seismic  disturbance  which  propagates  at  the  same 
rate  as  the  imposed  pressure  field.  It  may  be  readily  verified  that  this 
disturbance  is  identical  in  form  to  that  created  by  a  static  load  p  (x)  imposed 
upon  the  surface  of  the  half  space  (cf,  Sneddon,  1951).  Furthermore  it  is 
important  to  note  that  in  this  particular  case  there  is  no  energy  transport  by 
the  normal  seismic  waves  (i.e.,  waves  traveling  at  speeds  a  or  0  or  the 
Rayleigh  wave  velocity) .  To  an  observer  moving  at  the  same  speed  as  the  load, 
the  disturbance  would  remain  fixed  and  would  be  identical  to  that  created  by 
a  static  load. 


2.2  DISCUSSION 


The  basic  question  to  be  answered  in  this  technical  note  is  whether  or  not 
the  seismic  disturbance  created  by  the  pressure  wave  can  contribute  signifi¬ 
cantly  to  the  seismic  background.  In  order  to  answer  this  question,  the 
vertical  and  horizontal  displacements  created  by  a  sinusoid  of  angular  frequency 
o)0  traveling  at  a  speed  c0  meters/sec  with  a  pressure  amplitude  of  P  microbars 
are  determined  below  for  various  half-space  models.  For  this  case 


P(w) 


/ 


,i(u>-w0)t  ^ 


(44) 


=  211  P  6  (uj-tjo0) 


Using  equation  44  in  equations  42  and  43  and  performing  the  indicated 
integration,  the  displacement  components  are  found  to  be 


U(x,z,t) 


eiwo  (t-2E) 

_ 2a 


(45) 


W(x,z,t) 


(46) 


Now  at  the  surface  the  displacements  become 


U(x,o,t)  =  ^coP  e  ia)o  (t  "  cn)  (47) 

2(A+u)  w 0 

x 

=  CqP  _  (X+2y)  eia)o  (t  -  cQ  ) 

2p  (A+y)  a)0 


W(x,o,t) 


(48) 


tquations  47  and  48  have  been  used  to  compute  the  amplitudes  of  the 
displacement  components  for  various  models  with  c0  =  5  JJ_  -  u  <  n_ 

<  .  10  °  50 

(20  -  T  ■  100  sec)  and  P  =  i.  The  models  are  described  in  table  1.  The 

results  are  displayed  in  figures  1  through  4  for  both  the  vertical  and  hori¬ 
zontal  components  of  displacement. 


Table  1.  Parameters  of  the  models  used  to  compute  vertical 
_ and  horizontal  components  of  displacement _ 


a 

8 

P 

km 

km 

g 

Model  No. 

sec 

sec 

3 

cmJ 

Geology 

1 

0.3 

0.1 

1.6 

Unconsolidated 

sediments 

2a 

2.9 

0.7 

2.3 

Sedimentary  rock 

2b 

2.9 

1.0 

2.3 

Sedimentary  rock 

3 

4.8 

2.8 

2 . 65 

Sedimentary  rock 

4 

5.8 

3.25 

2.85 

Granites 

Model  1  is  clearly  an  extreme  case  and  was  included  to  show  the  worst  possible 
conditions.  As  shown  in  figure  1  this  model  is  quite  sensitive  to  small  changes 
in  atmospheric  pressure.  For  instance,  in  the  period  range  20-40  seconds,  the 
average  vertical  displacement  is  80  mu  for  a  pressure  amplitude  of  1  u  bar. 

In  the  same  range  the  average  horizontal  displacement  is  10  mp/p  bar.  These 
results  imply  that  one  of  the  reasons  that  long-period  installations  located 
on  thick  sections  of  alluvial  fill  are  generally  characterized  by  high  noise 
levels  is  that  the  unit  is  quite  sensitive  to  small  changes  in  atmospheric 
pressure . 

Models  2a-2b  have  velocities  and  densities  characteristic  of  many  sandstone 
shale  sequences  found  in  the  mid-continent  and  Gulf  Coast  regions  of  the 
United  States.  Model  2b  differs  from  2a  only  in  that  the  shear  wave  velocity 
has  been  increased  to  1  km/sec.  The  curves  shown  in  figure  2  indicate  that  the 
pressure  sensitivity  of  these  models  considerably  lower  than  the  uncon¬ 
solidated  sediments  model. 

In  the  period  range  20-40  seconds  the  average  vertical  displacement  is  about 
1  mp/p  ba*  for  Model  2a  and  0.6  mp/p  bar  for  Model  2b.  However,  if  we  assume 
pressure  amplitudes  on  the  order  of  100  p  bars,  the  corresponding  vertical 
displacements  of  the  surface  would  be  100  and  60  mp  for  Models  2a  and  2b.  Now, 
in  general,  the  vertical  seismic  background  observed  at  periods  greater  than 
25  seconds  is  somewhere  on  the  order  of  100-200  mp.  Thus,  the  pressure  waves 
with  amplitudes  of  100  p  bars  or  more  can  generate  p  significant  fraction  of 
the  seismic  background  at  periods  greater  than  25  seconds.  If  the  installations 
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AMPLITUDE 


PERIOD  (seconds) 


Figure  1.  Surface  displacement  response  of  a  half  space 
consisting  of  unconsolidated  sediments 
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AMPLITUDE 


1  1 

MODEL  NO.  2a 

I  1  1 

MODEL  NO.  2b 

a  =  2.9  X  103  METERS/SEC 

a  =  2.9  X  103  METERS/SEC 

.3  =  .7  X  103  METERS/SEC 

/?  =  i  n  y  m3  MFTPRR/^Fr 

p  =  2.3  g/cm3 

p  =  2.3  g/cm3  . 

C  =  5  METERS/SEC 

1  1 

C  =  5  METERS/SEC  / 

IU 


PERIOD  (seconds) 


Figure  2.  Surface  displacement  response  of  a  weak  sandstone 
shale  half  space 
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Figure  3.  Surface  displacement  response  for  a  competent 
sedimentary  rock  half  space 
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Figure  4.  Surface  displacement  response  for  a  granite  half  space 
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happen  to  be  located  on  a  thick  section  of  moderately  to  poorly  indurated 
sandstones  and  shales.  By  thick  it  is  meant  that  the  dimensions  of  section 
are  much  larger  than  the  wavelength  of  the  disturbance.  For  the  periods  and 
speeds  considered  above  the  corresponding  wavelengths  would  be  less  than 
200  meters.  Consequently,  these  results  could  be  applied  to  sections  whose 
dimensions  are  as  small  as  a  kilometer.  It  is  important  to  note  that  the 
shallow  structure  beneath  the  Montana  LASA  is  quite  similar  to  that  used  in 
Modei  2b  at  least  down  to  depths  of  about  1.5  km  (cf,  Glover  and  Alexander, 
19681.  Thus,  it  appears  reasonable  to  expect  that,  at  the  Montana  LASA, 
atmospherically-generated  seismic  noise  can  contribute  significantly  to  the 
vertical  long-period  noise  spectrum  as  suggested  by  Capon  (1968). 

The  situation  is  quite  different  for  the  mere  competent  rocks  which  are 
represented  by  Models  3  and  4.  As  shown  in  figures  3  and  4,  the  sensitivity 
of  these  models  to  atmospheric  pressure  fluctuations  has  been  reduced  by  about 
an  order  of  magnitude  relative  to  the  sandstone-shale  models  considered 
previously.  Therefore,  in  rocks  of  this  type,  displacements  caused  by  atmos¬ 
pheric  loading  cannot  contribute  significantly  to  the  long-period  seismic 
noise  background  unless  the  pressure  amplitudes  are  on  the  order  of  a  millibar 
or  more. 

The  previous  paragraphs  have  been  concerned  with  the  displacements  at  the 
surface  of  the  earth  caused  by  a  plane  pressure  wave.  These  waves  also  pro¬ 
duce  tilts,  which  can  be  a  major  source  of  noise  on  the  long-period  horizontal 
seismograph  system.  The  tilts  produced  by  a  plane  pressure  wave  may  be 
determined  in  the  following  way. 


3  (x,  z,  t)  be  the  angle  of  tilt  measured  from  the  horizontal.  Then 

dW 
dx 


Tan  0  (x,  z,  t) 


or  using  equation  46 

iP 


A+2. 


Tan  e  <x.z.t)  -  f  ei«o  (t  -  % 


(49) 


Since  0  is  generally  small 
Tan  0  %  0 

Consequently  at  z  =  0 


0  «  Tan  0 


P  A+2 
ZF  (A+y) 


ei  [wo  (t-  *  )  + 

co 


(50) 


Thus,  at  the  surface  the  amplitude  of  the  tilt  is  dependent  only  upon  the 
elastic  constants  of  the  medium  and  the  pressure  amplitude  and  is  independent 
of  frequency.  The  tilts  associated  with  a  1  microbar  pressure  amplitude  in 
each  of  the  four  models  described  above  are  listed  in  table  2. 
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It  can  be  seen  from  this  table  the  tilts  at  the  surface  are  quite  small. 
However,  because  of  the  tilt  sensitivity  of  the  long-period  horizontal  seis¬ 
mometer,  even  these  small  tilts  can  produce  large  apparent  horizontal  dis- 
displacements .  Rogers  (1968)  has  shown  that  the  apparent  horizontal  dis¬ 
placements  caused  by  tilt  is  given  by 

U°  =  0  417  G  (T)  (51) 


where  T  i s  the  period 

G  (T)  is  the  displacement  response  of  the  horizontal  seismograph 
system 

g  is  the  local  value  of  gravity 

Ug  is  the  apparent  trace  displacement  caused  by  tilt 

U,  was  computed  using  the  values  of  0  listed  in  table  2,  a  pressure  amplitude 
o£  100  microbars,  and  the  standard  LRSM  long-period  displacement  response  for 
20  *  x  -  100  seconds.  The  resulting  values  are  for  a  system  operating  at  a 
magnification  of  104  at  25  seconds  and  are  plotted  as  a  function  of  period  in 
figure  5.  The  curves  in  this  figure  indicate  that  during  periods  of  moderate 
atmospheric  turbulence  systems  located  in  thick  sections  of  unconsolidated 
sediments  will  become  virtually  unusuable.  The  situation  is  not  much  better 
for  the  case  of  Model  2.  In  a  thick  sandstone  shale  sequence,  the  tilts  would 
generate  trace  displacements  between  about  40  and  90  mm,  if  the  system  were 
operating  at  magnification  of  104. 


In  the  more  competent  rocks  a  pressure  amplitude  of  100  microbars  would  generate 
trace  amplitudes  varying  from  2-7  mm  of  trace  displacement  at  a  magnification 
of  104.  Thus,  it  would  appear  that  the  combination  of  moderate  atmospheric 
turbulence  and  the  tilt  sensitivity  of  the  long-period  horizontal  seismometer 
can  limit  the  magnifications  possible  at  the  surface  to  about  1  to  2  x  104 
even  for  the  more  competent  rocks.  This  is  in  qualitative  agreement  with  our 


experience  that  long-period  horizontal 
operated  at  magnifications  exceeding  2 


seismograph  systems  can  rarely  be 
x  104  at  or  near  the  surface  of  the 


earth . 
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TRACE  AMPLITUDE  FOR  A  MAGNIFICATION  OF  104  AT  25.6  SEC  (mm) 


0  20  40  60  80  100  120  140 


T  (seconds) 


Figure  5. 


Trace  displacements  on  a  horizontal  long-period  seismograph 
with  a  standard  LRSM  response  caused  by  earth  tilts 

G  5832 
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2.3  ATTENUATION  WITH  DEPTH 


An  interesting  feature  of  the  half-space  response  to  slowly  moving  plane 
pressure  waves  is  the  rate  at  which  the  disturbance  decays  away  with  depth. 
Let  Rw  (z)  denote  the  ratio  of  the  vertical  displacement  at  depth  to  the 
vertical  displacement  at  the  surface.  Then  from  equations  46  and  48 


C  A  2  u  +  z  .  u)0  - 

R  (z)  =  A+p  c0  j  e  c — 

w  - -  uo 

A+2p 

A+p 


(52) 


Similarly,  if  R^z)  denotes  the  ratio  of  the  horizontal  displacement  at  depth 
z  to  the  horizontal  displacement  at  the  surface,  then  from  equations  45  and  47 


V  +  Jfoz  u>o 

Hyd)  =  C0  «"  C0 

u 

A+p 


(53) 


If  we  let  R0  (z)  define  a  similar  ratio  for  the  tilts  it  is  clear  from  equations 
48  and  50  that 

Rg  (z)  =  Rn  (z)  (54) 


Rw(z)  is  plotted  as  a  function  of  z  in  figure  6  for  the  cases  where  c  =  5 
A  +  2p/A  +  p  =  1.0  and  1.5,  and  for  periods  20,  40,  and  60  seconds.  Generally, 
the  results  show  that  the  attenuation  with  depth  is  weikly  dependent  on  the 
elastic  properties  of  the  medium.  The  values  1.0  and  1.5  span  the  practical 
range  of  values  that  the  ratio  A  +  2p/A  +  p  can  assume  for  rocks,  with  the 
higher  values  being  associated  with  the  more  competent  rocks.  It  can  be  seen 
from  the  curves  in  this  figure  that  the  differences  in  attenuation  with  depth 
between  poorly  consolidated  and  competent  rocks  is  only  2  dB. 

Furthermore,  the  seismic  disturbance  will  decay  rapidly  with  increasing  depth. 
For  example,  at  a  depth  of  100  meters  (328  ft),  the  amplitude  of  a  20-second 
period  sinusoid  would  be  40  dB  below  the  amplitude  observed  at  the  surface. 

At  the  same  depth,  waves  with  periods  of  40  and  60  seconds  would  be  reduced  by 
about  11  and  6  dB,  respectively.  However,  in  order  to  obtain  a  40  dB  reduction 
of  the  noise  generated  by  atmospheric  pressure  variations  over  the  period 
range  from  20  to  60  seconds,  the  detectors  would  have  to  be  placed  at  least 
300  meters  deep. 
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Figure  6.  Attenuation  of  vertical  displacements  and  tilts  with  depth 


Thes*?  results  suggest  that  m  regions  where  local  atmospherically-generated 
seismic  noise  can  be  anticipated  to  be  significant  at  the  surface  (i.  e.,  in 
regions  of  thick  sandstone-shale  sequences  or  alluvial  fill)  much  of  this 

type  of  noise  can  be  eliminated  by  placing  the  detector  at  moderately  shallow 
depths. 

One  further  point  requires  emphasis.  Recently,  there  has  been  a  considerable 
interest  in  extending  the  response  of  long-period  systems  to  cover  a  pass 
band  from  30  to  100  seconds.  From  equations  47  and  48  it  is  seen  that  the 
displacement  response  of  the  half  space  increases  linearly  with  increasing 
period.  Furthermore,  equation  51  shows  chat  the  tilt  response  of  the  hori¬ 
zontal  seismograph  system  will  increase  as  the  square  of  the  period. 
Consequently,  on  the  basis  of  these  results,  one  can  anticipate  that  atmos- 
pherica„ly  generated  seismic  noise  will  place  severe  restrictions  on  the 
magnifications  that  can  be  obtained  with  ultra  long-period  systems  which  are 
operated  near  the  surface. 


3.  CONCLUSIONS 


The  basic  purpose  of  this  initial  study  was  to  determine  if  the  displacements 
and  tilts  generated  by  plane  pressure  waves  can  contribute  significantly  to 
the  long-period  noise  field.  On  the  basis  of  the  results  presented  in  the 
previous  paragraphs,  it  appears  that  the  vertical  displacements  generated  by 
these  waves  can  contribute  significantly  to  the  long-period  noise  field  pro¬ 
vided  that 

a.  The  installation  is  located  at  or  near  the  surface  of  a  section  of 
alluvial  fill  or  poorly  to  moderately  indurated  sandstones  and  shales  which 
is  at  least  1  km  thick,  and 

b.  The  pressure  amplitudes  are  of  the  order  of  100  microbars  or  greater. 

Pressure  waves  with  amplitudes  of  the  order  of  100  microbars  will  also  create 
significant  tilt  noise  on  long-period  horizontal  seismograph  systems  located 
at  or  near  the  surface,  regardless  of  the  rock  type.  However,  the  noise  will 
be  much  less  severe  on  competent  units. 

The  results  also  show  that  the  seismic  disturbances  created  by  plane  pressure 
waves  will  decay  rapidly  with  depth.  Furthermore,  the  rate  of  decay  is 
virtually  independent  of  the  rock  type.  The  implications  are  that  the  contri¬ 
bution  of  atmospherically-generated  noise  to  the  long-period  background  can 
be  virtually  eliminated  for  periods  less  than  a  minute  by  placing  the  detectors 
at  depths  on  the  order  of  300  meters. 

Finally,  it  appears  that  atmospherically  generated  noise  will  be  more  severe 
on  ultra  long -period  systems  than  on  presently  operating  systems. 

These  results  are  considered  to  be  sufficiently  significant  to  warrant  the 
study  of  the  response  of  more  complex  media  to  more  realistic  pressure  fields. 
Consequently,  numerical  techniques  to  compute  the  response  of  layered  media 
to  arbitrary  pressure  fields  are  being  developed. 
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STATIC  LOADING  OF  A  LAYERED  ELASTIC  HALF  SPACE  BY  ATMOSPHERIC  PRESSURE 
_  _  VARIATIONS 


1.  INTRODUCTION 


The  response  of  a  homogeneous  elastic  half  space  to  moving  loads  has  been 
discussed  in  Appendix  1  (Sorrells) .  This  appendix  shows  how  to  compute 
the  response  of  a  layered  elastic  half  space  to  static  loads.  The  Haskell- 
Thompson  matrix  method  is  used  to  derive  displacements  and  tilts  for  the 
two  dimensional  case.  It  is  likely  that  the  layered  half  space  is  a  more 
satisfactory  model  in  actual  field  situations,  and  the  computer  programs 
will  be  used  to  verify  the  theory  in  conjunction  with  field  experiments. 


2.  GENERAL 

S  ppose  that  the  stresses  at  the  free  surf,  ce  of  a  layered  elastic  solid 
(the  earth)  due  to  atmospheric  pressure  variations  can  be  described  as 
space- time  stationary  random  processes  with  power  spectral  densities  P3i(w>  kj,  ^2) 
where  w  is  the  angular  frequency  and  k^  and  ky  are  the  horizontal  wave 
numbers  and  i  =  1,  2,  3  indices  for  the  three  coordinate  axes.  It  is 
assumed  that  the  coordinate  axes  Xj  and  X2  are  horizontal  and  X3  is  pointing 
downward  into  the  elastic  half  space.  If  the  velocities  of  the  pressure  vari¬ 
ation  are  considerably  below  the  elastic  velocities  in  the  elastic  half  space, 
the  deformations  of  the  half  space  can  t?e  considered  as  static  and  the  depend¬ 
ence  of  the  response  on  u  can  be  neglected.  Therefore,  the  response  of  the 
elastic  half  space  will  depend,  in  the  first  approximation,  only  on  the  hori¬ 
zontal  wave  numbers  k^  and  k2-  Both  shear  and  normal  stresses  can  act  on  the 
surface.  Shear  stresses  are  usually  due  to  wind  drag  and  pressure  on  irregu¬ 
larities  ard  obstacles  at  the  surface,  such  as  trees,  hills,  buildings,  etc. 

In  the  following  discussion  it  will  be  assumed  that  the  earth's  surface  is 
flat,  eJid  the  seismometers  are  located  far  from  obstacles  which  would  cause 
wind  drag.  For  a  flat  earth  the  normal  stresses  P33  (w,  kx,  k2)  would  cause  the 
tilts  observed.  In  order  to  derive  the  pressure  variations  observed  by  a 
microbarograph  it  is  assumed  that  the  components  of  the  pressure  field  with 
different  kj  and  k2  are  independent  and  therefore  the  power  spectra  of  these 
pressures  are  additive.  This  is  not  an  unreasonable  assumption  and  therefore 
power  spectra  of  the  microbarograph  output  can  be  obtained  by  integration  over 
ki  and  k2  giving 


^mm  0*0 


P33(w,ki,k2)dk1dk2 


—00 


(1) 
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wh.-re  , 1(a)  is  the  frequency  response  of  the  microbarograph.  A  horizontal 
seismograpn  oriented  in  a  direction  Xj  located  at  a  certain  depth  d,  below 

logic"  ^  thC  Ult  With  the  power  sPectra  Cus ing  the  same 

X 

jPss(u)  s  ^  f  f  Ij2^kl>k2>dl)p33^a)»kl»k2)dklcik2  (2) 

— x 

Where  S(m)  is  the  seismometer  tilt  response  function  and  Ti (k, ,ki,d,)  is  the 

*‘t  cr  l“  ,reSp°nse  of  the  sJound  at  dePth  dl  t0  normal  ljads  of  the  form 
exp  i  (h,x,*k2x2)  acting  at  the  surface.  The  cross  spectra  between  the  m'rro- 

a^Ss  d  ™dT0ntual  seif”eter  or  two  horizontal  seismometers  located 
at  depths  d^  and  d2  and  similarly  oriented  will  be 

ipsm(w)  =  M(a))S(a))J^y  Tj(k1,k2,d1)  P33 (oi ,kj , k2)  dkjdk2  (3) 

jPs^^)  =  S1(w)S2(a) )  f  j  Tj(k1>k2,d1)Tj(k1,k2,d2)P33(a),k1)k2)dk1dk2 

where  S^w)  and  S2 (w)  are  the  responses  of  seismometers.  Similar  expressions 
can  be  written  for  the  response  of  the  vertical  seismometers,  but  since  these 
instruments  are  sensitive  to  vertical  displacements  and  not  to  tilts  the 
effect  of  atmospheric  loading  would  be  much  smaller  on  these  instruments  It 

can  be  shown  from  the  above  expression  that  even  if  the  noise  from  other 
sources  is  absent,  and  the  atmospheric  noise  is  propagating  unidirectional ly 
the  coherence  between  the  outputs  of  the  microbarograph  and  a  seismometer  or 
between  two  seismometers  located  at  two  different  depths  is  always  less  the 

haifsDace'at  ~  C°n!t’  ^  tiU  response  is  Constant  for  a  homogeneous 

halfspace  at  x3=0  for  the  two-dimensional  case  where  one  of  the  horizontal 

wavenumbers  is  zero.  If  such  model  of  the  subsurface  and  unidirectional 

propagating  noise  is  applicable,  the  coherence  between  the  seismometer  tilt 

output  and  the  microbarograph  may  be  close  to  unity.  In  other  cases  it  is 

variations t0  ^  f’R  structure  of  the  atmospheric  pressure 

In  order  to  determine  the  f-k  structure  of  the  atmospheric  pressure  variations 
microbarograph  array  measurements  are  necessary.  If  such  measurements  are 
available  the  above  theory  can  be  verified  from  simultaneous  pressure  and 
tilt  recordings.  r-(k1,k2,x3)  can  be  determined  numerically  if  the  elastic 
parameters  of  the  subsurface  are  known.  It  is  assumed  here  that  the  subsurface 
can  be  modeled  with  a  layered  elastic  half  space,  and  the  effect  of  the  vault 
cavity  at  depth  can  be  neglected. 

Ihe  purpose  of  this  study  is  to  develop  methods  to  eliminate  tilt  noise  caused 
by  atmospheric  pressure  variations  on  horizontal  long-period  seismometers  If 
the  f-k  structure  of  the  pressure  variations  is  such  that  a  reasonable  number 
of  microbarographs  is  sufficient  to  predict  the  tilt  noise  on  the  seismometers 
a  significant  reduction  of  the  noise  on  long-period  seismographs  is  possible. 

A  single  microbarograph  is  not  likely  to  achieve  this  goal,  since  the  f-k  field 
epends  on  the  wincj  direction  and  even  if  a  halfspace  model  and  two-dimensional 
pressure  field  are  applicable  the  sign  of  the  phase  lag  of  ir/2  between  the  pres¬ 
sure  variations  and  the  tilt,  and  the  magnitude  of  tilt  in  direction  xi  must  be 
determined  from  the  horizontal  pressure  gradient.  To  estimate  this  gradient  an 
array  of  at  least  three  microbarographs  is  necessary.  The  knowledge  of  wind 
direct .on  and  speed  may  reduce  the  number  of  microbarographs  necessary. 
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If  the  average  f-k  structure  of  the  atmospheric  pressure  variations  is  known 
at  a  certain  location  one  can  estimate  the  number  of  microbarographs  and 
seismometers  needed  to  eliminate  a  certain  fraction  $  of  loading  noise  power 
by  computing  the  multiple  coherence  functions  y(w),  defined  by  (Bendat  and 
Piersol,  1966)  as 

Y2M  -  1  -  [GiMGV)]'1, 

for  various  combinations  of  sensors,  where  Gi(w)  and  G1  (to)  are  the  i-th 
diagonal  elements  of  the  spectral  matrix  and  inverse  spectral  matrix  of 
sensor  outputs.  The  elements  of  these  matrices  can  be  computed  for  any 
layered  structure  and  instrument  combination  by  using  formulas  of  the  type 
(1),  (2),  and  (3).  The  sensor  combination  which  yields  Yj(uO  >_  B  would  be 
sufficient  to  eliminate  this  fraction  of  loading  noise  power  at  the  i-th  sensor. 

The  mathematics  used  in  computation  of  tilt  responses  is  described  in 
appendices  2a  and  b.  The  computer  programs  are  described  in  appendix  2c. 
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PROPAGATION  MATRIX  FOR  TWO-DIMENSIONAL  DEFORMATION  OF  A  LAYERED  ELASTIC 

_ half  space  due  to  surface  line  stresses 


This  theory  is  based  on  the  extension  of  results  described  by  Sneddon  (1951) 
for  two-dimensiona1  deformation  of  elastic  solids,  the  Fourier  transform  G 

°  stress  Action  with  respect  to  the  horizontal  coordinate  x,  is 

satisfied  by  the  equation  1 


(1) 


where  k  is  the  horizontal  wave  number,  x,  is  the  vertical  coordinate  This 
equation  has  the  general  solution 


G  =  (A  +  Bx.j)e 


(C  +  Dx3)e 


(2) 


The  Fourier  transforms  of  the  stresses  and  displacements  can  be  derived  from 
this  equation  the  following  way 


(3) 


wh®re,  u»  v»  Txx>  Txy  are  the  vertical  and  horizontal  displacements,  normal 
and  sheer  stresses  at  the  surface,  respectively,  0  is  the  Poisson  ratio,  E  is 
the  Young’s  modulus  of  the  elastic  medium. 

For  a  layered  medium  consisting  of  homogenous  isotropic  layers  the  Solution 
(2)  can  be  propagated  by  the  Thompson -Haskell  matrix  method  where  ABC 
and  D  take  the  place  of  constants  A’  ,  A’’m,  w’m,  A"m  in  Haskell’s  paper!  The 
matrix  which  is  equivalent  of  Haskell’s  Dm  for  this  boundary  value  problem 
can  be  obtained  by  substituting  2  into  3,  It  takes  the  form 
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INTEGRATION  IN  THE  kx  -  k2  PLANE 


INTEGRATION  IN  THE  kx  -  k2  PLANE 


The  two-dimensional  loading  program  can  be  used  to  compute  the  integrals  of 
the  type  (1),  (2),  and  (3)  if  the  integration  performed  is  a  polar  coordinate 
system.  In  this  case  k^  =  k  cos$,  k2  -  k  sind,  dk^dk2  =  kdkdcf  ~ 

P(«,k,$),  T1(k1,k2,d)  -  T(k,d)  cos<f>,  etc. 


Therefore,  the  power  spectrum  of  the  microbarograph  will  become 

„  2n 


pnunW  *  «'  W 


// 

k=0  <j)=0 


k  P 33  > k  ,  4> )  d4>dk 


The  tilt  power  spectrum  on  the  seismometers  oriented  in  the  directions  of 
horizontal  coordinate  axis  and  X2,  respectively,  and  located  at  depth  d 
will  become 


Pll(u)  =  Si(u>) 


If 

k=0  4>=0 


k  P33(uj,k,b)T^(k,d)cos^4i  d4>dk 


»  2n 

P22  (u))  =  S2^  /  /  k  P33^»k>^  T2(k,d) sin24>  d^dk 

k=0  4)=0 


where  <J>  is  measured  from  the  x^  axis 

The  cross  spectrum  between  the  microbarograph  and  the  seismometer  oriented  in 
the  Xj  direction  will  become  2ji 

Pml  ^  =  Sl^  f  j  k  P33^U)#k#(J)^  T(k'd)  lCOS(H  d(frdk 

k=0  d=0 


A  number  of  similar  formulas  can  be  derived  for  other  instrument  combinations. 


-1- 


TR  70-12,  app  2b 
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GENERAL  DESCRIPTION  OF  COMPUTER  PROGRAM  LOAD  2 


i 


- 


GENERAL  DESCRIPTION  OF  COMPUTER  PROGRAM  LOAD  2 


A  computer  program  LOAD  2  has  been  written  for  the  calculation  of  the 
displacements  and  tilts.  The  program  computes  the  responses  of  the  layered 
medium  to  sinusoidal  stresses  at  the  surface  as  functions  of  wavelength  and 
depth.  Both  normal  and  shear  stresses  can  be  applied.  The  program  reads 
the  densities,  the  stresses  and  seismic  velocities  of  the  layers,  the  ranges 
of  horizontal  wave  numbers  at  the  surface.  It  is  a  simple  matter  to  use  the 
program  to  compute  the  deformation  of  the  layered  medium  due  to  stresses  of 
arbitrary  shape  by  Fourier  synthesis  of  sinusoidal  loads.  The  same  program 
can  be  used  to  compute  strains  due  to  surface  stresses  at  any  depth,  and  with 
some  modifications,  to  compute  stresses  and  displacements  due  to  dislocations 
of  certain  types  in  a  layered  medium. 

Some  numerical  problems  can  arise  if  the  displacements  are  computed  for  greater 
depths  using  large  wave  numbers;  after  reading  the  value  near  zero  at  a  certain 
depth  the  displacements  start  to  increase  without  bound  with  increasing  depth. 
The  problem  is  purely  numerical  and  is  similar  in  nature  to  that  encountered 
in  surface  wave  calculations  (Dunkin  1965) .  It  can  be  remedied  by  double  pre¬ 
cision  calculations  or  the  elimination  of  the  large  (and  obviously  erroneous) 
displacement  values  below  a  certain  depth. 

The  matrix  alegebra  is  analogous  to  that  presented  by  Haskell  (1953),  products 
like  Dn-Ejj,-1  are  computed  numerically  by  inverting  Em  and  multiplying  instead 
of  deriving  all  the  algebra. 

The  program  has  been  checked  out  against  formulas  derived  by  Sorrells  (1969) 
and  several  other  examples  in  the  literature. 

Figures  1  to  6  show  some  examples  of  calculations  performed  using  the  formulas 
presented  above. 
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Figure  3.  Half-space,  depth  of  25  meters 
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Ii8fkeS  k  t0H?  Sh°f  C“t0Ur  PlT  0f  the  Ult  r«P°nses  in  the  x,  direction 
face  1  d>  functl°[ls  of  horizontal  wave  numbers  kx  and  k2  for  a  sur- 

used  Li  Life™  eXP  (  1  "  1  *  k2  x  2)  with  1  ubar  amplitude.  The  model 

km/sec.  respectiveIy!1VZ5rity120n51g/»3She?hreVk11itre5  °f  S'5  “d  3’3 
for  three  differentlepths;  0  (slrfacef^u'andlsIi^ThlloTtLI  arl” 

,  ,  d  intervals  and  the  tilt  is  given  in  radians.  The  contours  were  done  by 

a  uraler  hHlh  "f  °f  n“erical  val“«  and  therefore  are  not  very 
accurate,  but  the  figures  are  given  only  as  illustrations  of  the  general 

nature  of  the  responses.  On  the  k  axis  the  response  has  to  be  zero  Ihe 
figures  show  the  fast  decay  of  the  short  wavelengths  with  depth.  Figures  4  to 
show  the  tut  responses  for  a  model  with  a  low  velocity  layer  overlying  a 
half  space  with  elastic  constants  given  above.  y  8 

The^velocities  and  density  for  thb  low  velocity  layer  are  the  following- 

liver  islsl’  V'J'6  k1Sec  and  6  =  2‘  g/™3-  The  thickness  of  thf’ 

th„r  f  2u  Ti  A  the  surface  for  large  wave  numbers  the  response  resembles 

layerl  ^  Wlth  elastlc  constants  equalling  those  of  the  low-velocity 

we?«h  the  tUt  iS  very  snla11  and  11  is  almost  equal  to  the  til<-  response 
for  the  half  space.  This  indicates  that  burial  of  seismometers  is  a  very 

1  techna<>ua  for  teducing  tilt  noise.  If  the  seismometer  is  buried  close 

a?  ™!r/' "1“-  0f  3  aurface  low-velocity  layer  the  tilt  noise  is  about 

at  the  same  level  as  it  would  be  in  competent  rocks  at  the  same  depth. 
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IJ  ABSTRACT 

^Preliminary  studies  of  the  long-period  .seismic  noise  and  the  atmospheric  pressure 
fields  were  made  at  the  Kleer  mine  near  Grand  Saline,  Texas,  during  the  winter  of 
1968.  The  results  of  the  studies  indicate  that  during  periods  of  atmospheric 
turbulence  the  vertical  and  horizontal  components  of  the  seismic  noise  field 
observed  at  the  surface  may  increase  in  total  power  by  as  much  as  16  and  $4  dB 
within  the  period  range  from  20  to  60  seconds.  In  contrast,  long-period  systems 
operated  in  the  mine  show  little  change  in  power  levels  regardless  of  the 
atmospheric  conditions  at  the  surface.  The  observations  are  consistent  with  the 
hypothesis  that  the  increase  observed  at  the  surface  is  caused  by  atmospherically- 
generated  deformations  of  the  earth.  However,  because  of  questions  regarding  the 
surface  vaults  capability  to  attenuate  pressure  changes,  the  possibility  that  the 
noise  was  caused  by  vault  leakage  must  be  considered  an  equally  likely  explanation 
at  the  present  time.  The  results  indicate  that  atmospherically-generated  seismic 
noise  can  seriously  degrade  the  capability  of  a  long-period  surface  installation 
to  detect  surface  waves.  Because  of  this  finding,  further  more  definitive  studies 
are  recommended. 
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